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Abstract 

We study a novel class of nonsingular time-symmetric cosmological 
bounces. In this class of four dimensional models the bounce is induced 
by a perfect fluid with a negative energy density. Metric perturbations 
are solved in an analytic way all through the bounce. The conditions for 
generating a scale invariant spectrum of tensor and scalar metric pertur- 
bations are discussed. 



1 Introduction 

Bouncing universes become attractive and fashionable periodically. Bouncing 
models can solve the horizon problem with a long period of slow contraction 
(a ~ {—tY with < p < 1) as opposed to inflationary models which undergo a 
period of quasi-exponential expansion. Models with a bounce play an important 
role in string cosmology (see [1] and references therein) . It is therefore interesting 
to study how metric perturbations are generated in bouncing models paying 
particular attention to what happens close to the bounce, even in a simplified 
model. 

In this paper we study a hydrodynamical toy model where the bounce is 
induced by a negative energy perfect fluid. Negative energy density fluids seem 
ideas quite wild - or desperate [2] , but they have already been considered as toy 
models [3, 4]. Moreover, negative energy fluid have been recently connected to 
brane models [5], where the interaction between bulk and brane can induce a 
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negative energy component in the brane through the projection of the 5d Weyl 
tensor. 

A two field description of a bounce is also very pedagogical for the treat- 
ment of cosmological perturbations through a bounce. In such a bounce, both 
the Hubble parameter H and its time derivative H cross zero. Only the latter 
phenomena occurs in an expanding universe. Indeed, the Hubble parameter 
changes slope during the transition from an accelerated to a decelerated expan- 
sion (from inflation to the standard big bang era, i. e. during reheating) and 
viceversa (from the matter dominated era to the present accellerated stage). 
The vanishing of H is instead a peculiarity of a bounce. It turns out that the 
treatment of cosmological perturbations during a bounce is much more compli- 
cated than during a stage which only includes H = 0. 

The paper is organized as follows. In section 2 we introduce a novel class of 
bouncing universe. To our knowledge this class generalizes the case of a radiation 
dominated universe where the bounce is obtained with a massless scalar field 
with negative energy density [2, 3, 4]. In section 3 we study both in analytic 
and numerical ways the evolution of gravitational waves in such a background. 
The equation can be solved analytically with the use of spheroidal functions, as 
Sahni did for gravitational waves in a radiation-matter transition [6] . In section 
4 we describe the problems in treating cosmological scalar perturbations during 
a bounce. In section 5 we study the evolution of scalar fluctuations and in 
section 6 we conclude. 



2 Background 

We consider a bouncing universe which satisfies the Einstein equation: 
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where SttG = -Mpi^, m < n and p+, p_ are constants. Prom this last relation it 
is obvious that p_ is important only close to the bounce. The two fluids have 
ratios between pressure and energy densities w+ = m/3 — 1 and W- = n/3 — 1, 
respectively. 
By assuming 

n = 2(m - 1) ^ u;_ = 2w+ + ^ (2) 
the scale factor a in terms of the conformal time is: 
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This class of solutions covers many different cases of physical contraction, 
each of them bounced in an expansion phase by a negative energy density fluid 
with a different equation of state. The parameter e regulates the magnitude of 
the scale factor at the bounce and the time scale r]o sets the duration of the 
bounce. Contraction driven by dust, radiation, free scalar field and ultra-stiff 
equation of state {w+ >> 1) arc given by q: = 1,1/2, 1/4 and a ^ 0, respectively. 
Since W- > the bounce is unavoidable in Eq. (1) and it does not depend on 
initial conditions. 

Important relations are: 
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where x = rj/rjo- We observe that p- cannot be simply related to the cur- 
vature of the spatial sections since m > 3 to avoid negative instability for fluid 
fluctuations ^ . However, the above solution can be extendend to non flat spatial 
sections, as for the particular case a = 1/2 [5]. The Hubble law in Eq. (1) in 
general is: 

1 r „ ^ -, 

(11) 
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where K = 0,1, —1 corresponds to flat, closed and open spatial sections. For 
K = 1 the solution for the scale factor is: 



a{r]) 



P+ 
6M2 



1 — ci cos 



(12) 



^Scalar metric fluctuations have been studied during a bounce in [10] for a universe filled 
by a massive scalar field with closed spatial sections. 
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Figure 2: The curvature term a"/a as a function of a; for a = 1/2,1 and rio = 1. 
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Figure 3: The curvature term a!' ja as a function of x for a = 0.1,0.2 and 




Figure 4: The curvature term a" ja as a function of a and x. 
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where cx = Jl — 12Mpjp_/p^. For K = —1 the solution for the scale factor is: 



a(ry) 



P+ 



C2 cos h 



(13) 



where C2 — y 1 + 12Mpjp_//c)^. Also for the open and closed solutions a is 
given by Eq. (5). 



3 Gravitational Waves 

As usual, the amplitude of gravitational waves h satisfies the following equation: 

h" + 2nh' + k'^h = 0. (14) 
By neglecting the infinite wavelength solution is: 

h^A + B J ^=A+^Bx2F,(^^,2a;^;-x'^ (15) 

In particular, for a dust dominated collapse (a = 1) we get 

770 / X arctana;\ 

"mat ~ ^mat + -«mat ^ \^2(a;2 + 1) 2 ) ' ^ ' 

and for a radiation dominated collapse (a = 1/2) 

hrnd ~ ^rad + -Brad ^ arctana; , (17) 
Instead, for a contraction driven by a free scalar field, we get: 

hff, 1^ + arcsinha; . (18) 

The presence of at the denominator indicates that the amplitude of grav- 
itational waves blows up in the limit of singular bounces, i. e. e 0. This 
fact is not surprising and is related to the normalization condition between the 
two independent solutions for the modes W{h, h*) = i/a^. Fluctuations become 
singular in a singularity. 

It is interesting to note that the B mode (odd with respect to a; ^ —x) 
tends to a constant far from the bounce, only for a > 1/2. For a < 1/2 
the amplitude of gravitational waves grows outside the Hubble radius, even 
in the expanding phase. This result is not completely surprising, since for 
slow expansion {a{r]) oc {rjy with p < 1/2), free scalar fields (and henceforth 
gravitational waves) are a superposition of a constant and growing mode (but 
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Figure 5: Behaviour of a long-wavelength gravitational wave across the bounce. 
The plots are for a = .1,. 25, .5.1 (from bottom to top). 



proportional to krj). Here one has this effects for a different range of parameters, 
probably because of the bounce. 

These results indicate that it is potentially important to know the k depen- 
dence of B, in order to predict the spectrum of gravitational waves far from the 
bounce. The contribution of the B mode to h right at the bounce is zero in this 
infinite wavelength approximation. 

With the evolution of the scale factor in Eq. (3), the equation for h can be 
solved analytically. The trick is to put the above equation in the form of the 
differential equation of spheroidal wave function [7]. It is useful to introduce 
then: 

x = r]/T]o h = h{l + x^)"-^^^ . (19) 
After a little algebra one gets: 



dx'^ 



2x dh 
1 + x^ dx 
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We note that the above equation for krjo — reduces to the Legendre differential 
equation [7, 8] . The solution to Eq. (20) is given by the radial oblate spheroidal 
function [8]. The above equation can be obtained by the transformation z — > 
±ix in the equation in prolate coordinates [7, 8, 6]: 
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(21) 



where jj, ,6 are numerical coefhcients. In many books only the case with 
integer /i , is treated [7, 8], however fi ,i> ,j can be even complex numbers [9]. 
The solution for h is: 



.y^(i + x^)V-«^-:J(3)(,,^), 



(22) 
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Figure 6: Behaviour of a the phase of a long- wavelength gravitational wave 
across the bounce for a = 1. 



where (the following series is convergent only when \z\ > 1 [9], to have conver- 
gence at all times one needs a different expression): 

7 = —ikr]o , z = —ix , 72; = —kr] (24) 



^ = 1/ = 2a - 1 (25) 

The normalization factor is fixed by requiring that for a; —> 00 (and then 

\kr]\ 00): 



a\/2k 

1/2 



(26) 



Jim5.^(3)(x,7)=(^) i?«/,(7.) (27) 

In order to understand the behaviour for large wavelengths we follow Sahni 
[6] and we take the limit for 7 ^ ^ [7] : 

lim (3) {x, 7) = 7"/' {*iPi:{x) + #2Q'i{x)) + ^Qiiix) (28) 

where [x) , Q'^ {x) are respectively the associate Legendre functions of the first 
and second kind [7] and #i,#2,#3 are numerical coefficients. Gravitational 
waves have finite amplitude right at the bounce since -P^(O), (3(^(0) are finite 

^This limit identifies the infrajred modes at the bounce, i. e. k'^ « a" /a at r; ~ 0. 
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Figure 7: A snapshot of the gravitational wave amphtude as a function of 
logiokrjo at a; = 30 for a= 1/2. The other detail of the simulation are Xi = —10^, 
e = 10^ and a initial vacuum spectrum for gravitational waves. The linear fit is 
-0.8 - 0.5(Zofifiofcr/o + 3.4). 



(see Eqs. 8.6.1 and 8.6.2 in [8]). This leads to a spectrum for gravitational 
waves in the infrared: 

lim ~ —1— . (29) 

The above result is in formal agreement with Sahni [6] . 

Wc note that the phases of gravitational waves arc not constant while outside 
the Hubble radius, as shown in Fig. (6). Phases are locked to a constant value 
when the wavelength of a mode exceeds the Hubble radius during inflation: this 
phenomenon has been dubbed "dccohcrcncc without dccohcrcncc" [11]. Models 
based on a contraction which aim to produce a scale invariant spectrum seem 
to not have this peculiarity. 

The final spectrum of gravitational waves far from the bounce is scale invari- 
ant for a dust contraction (a = 1), and it is a vacuum spectrum for a radiation 
contraction (a = 1/2), as we can see from Figs. (7,8). 



4 Equations for Scalar Perturbations 

We now review the equations of motion of scalar perturbations in the longitu- 
dinal gauge [12] with two hydrodynamical fluids. The energy and momentum 
constraints are respectively: 

-A<E> + m ($' + n^) = A-kGo^ {6p+ - dp-) (30) 

+ W$ = AttG— [p+(l + w+)v+ -p-(l + w-)v-] (31) 
k 

The spatial part of the Einstein equations lead to: 

+ 3W$' + [2H' + H^]^ = A-kGo^ {6p+ - Sp-) , (32) 
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Figure 8: A snapshot of the gravitational wave amphtudc as a funetion of 
logiokrjo at a; = 30 for a = 1. The other details the same of the previous figure. 
The linear fit is 2.78 - 1.5{logiokr]o + 3.4). 

Usually Eq. (30) times the total sound speed c^: 

= P = P+~P 
^~ P P+- P 

is subtracted by Eq. (32) to obtain: 

+ 3(1 + cl)H^' + [c^fc2 + 2H' + + 2,cln'^] $ = ^-kGo^SS , (34) 
where 6S are the total non adiabatic pressure perturbations: 

6S = J2^Pi-'^lJ2^Pi- (35) 

i i 

In problems with a bounce the sound speed as defined in Eq. (33) becomes 
singular ^. This is connected with the violation of the null energy condition 
(henceforth NEC) p + p>0, which occurs not right at the bounce (?7 = 0), but 
when 

^ = -2i^(^ + ^^ = -2]^I^(^ + "'^^' 

pi pi i 

vanishes, i.e. when = 770/2 in this class of models. 

This problem has been already encountered in the study of cosmological per- 
turbations during reheating [14]. In a scalar field driven cosmology H vanishes 

when the field bounces, i.e. (/) ~ 0, otherwise it remains always negative. The 
equation for $ has singularities [14]. In [14, 15] the problem was approached 

•^Recently a paper which studies numerically the spectrum of <E> fluctuations during a fictions 
bounce driven by a single perfect fluid appeared [13]. We think that the toy model used in 
[13] is inconsistent. Infact the equation of motion of <1> used in [13] contains a regular and ad 
hoc quantity which substitues the physical c^, which becomes singular during the bounce. 
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(33) 



by using the Mukhanov variable Q{= Scj) + (j)'^ / H) whose evolution is regular 
during the oscillation of a scalar field, despite the NEC violation. The regularity 
of the equations which involve the Q& persist also in the multifield case [15]. 

When a bounce is present the equations governing the evolutions of the Qs 
become singular in i? = (see [15] for scalar fields). This would be true also 
for perfect fluids. Infact, for a system of perfect fluids characterized by Wp and 
c|,, the Mukhanov variables Qps would satisfy: 



Qf + SHQf 



'^H{l-v:p) + ^H\l-w%) 



Qi 



^ X [Qs, Qs of other components) . (37) 
H 



We note that in presence of just one component, only a dust contraction can 
generate a scale invariant spectrum, as for the scalar field case [16]. 
We therefore turn to different variables in the next section. 



5 Evolution of Scalar Perturbations 

In the appendix of [4] Peter and Pinto-Neto suggest a possible way to deal with 
metric perturbations during a bounce. They suggest to split the Newtonian 
potential $ in two components 

$ = (38) 

where each one satisfies: 

+ 3(1 + wi)n'i>- + [wik'^ + 2n' + n'^ + swiU'^] $i = o (39) 

with i = + ,—. This splitting is valid under the assumption that non adiabatic 
pressure perturbations in the two components are absent: 

5pi = cfSpi with = Wi (40) 

It is important to note that the above equation(s) are regular and uncoupled. 

It is very interesting to see that the contribution of to $ decays in the 
expanding phase. The reason for this behaviour is that $_ docs not have a 
constant mode in the long wavelength limit far from the bounce, as instead 
does. In order to see this, it is sufficient to evaluate 

2n' + n^ + zwiH"^ = -z{w - Wi)n'^ . (41) 

For i = — we get: 
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Figure 9: Behaviour of a long-wavelength mode of $_ across the bounce. After 
a transient decay around the bounce, <&_ does not stabilize at a constant value, 
as also seen analytically. The plots are for a = 0.1,0.25,0.5, 1 (from bottom to 
top after the bounce). 



which is a strictly positive quantity at all times. For i = + we get a vanishing 
quantity at leading order far from the bounce: 

The above result implies that the growing mode of during the contracting 
phase does not match to the growing mode of <i> in the expanding phase. This 
result is confirmed by numerical analysis, as shown in Fig. (9). 

Therefore, the — component is important only close to the bounce both at 
the background level and at the linear level. 

We now turn to the analysys of and we repeat the previous treatment 
for long wavelengths. For i = + we get a vanishing frequency at leading order 
far from the bounce: 

Now we would like to recast Eq. (39) in a spheroidal equation form. By intro- 
ducing: 

u, = {l + x')-^^, (45) 
where 5 — 1/2 — 3a(l + Wi)/2. The equation for Ui is: 

/, 2x / 



' 1+X2 ' 

where 



7 2 2 



U^^Q. (46) 



l + x'^ (1+^2)2 _ 

6i = 4a(a(l -f 3w0 - 1) + 6a(5(l + Wi) (47) 
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Figure 10: Behaviour of a long-wavelength mode of $+ across the bounce. After 
a transient decay around the bounce, <i>+ stabilizes at a constant value, as also 
seen analytically. The plots are for a — 0.1,0.25,0.5.,! (from bottom to top 
after the bounce). 

c, = 4a(l + 25) + 4(5^ (48) 
When i = +, we get 6+ ~ —2a{2a + 1), c+ = 4a(l + a), leading to: 



iy+ ^2a fi+ = y/4:a{l + a) . (49) 

The solution is therefore 



'J'+ = (r^^-^'^(-'^) (50) 



where M ^ 0{'-f ^/^) is fixed by imposing: 



on short wavelengths far from the bounce. Indeed, 

1/2 



lim lim = lim 

'kri — >oo X — yoo —krj — ^oo 

N / 1 



6-^'=" (52) 



(1 + a;2)" 

For large wavelengths, by using Eq. (29) one has: 

(53) 

A nearly scale invariant spectrum for emerges in the expanding phase for a ^ 
0. For all the other relevant values of a the spectrum for in the expanding 
phase seems too red to agree with observations (for a = 1/2 ~ fc~^/^ in 
agreement with [4]). 
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6 Discussions and Conclusions 



We have studied cosmological perturbations through a non-singular bounce in 
a simple class of with two perfect fluid toy models. A normal fluid drives the 
contraction/expansion, while a negative energy density fluid models the bounce. 
The equations of state of the two fluids are linked in a particular way described 
by Eq. (2). 

The study of scalar metric fluctuations during a non-singular bounce is rather 
tricky, but it can be approached, even analytically, as this example shows. Fol- 
lowing [4] , one could find variables whose evolution equation is regular through 
the bounce. We do not think that the approach proposed in [4] is unique in 
order to evolve cosmological scalar perturbations through a bounce. One could 
define a fictions total sound speed which is regular through the bounce and 
asympotically behaves like the sound speed of the dominant component. There 
is a possibility that in this way also the equation of motion for isocurvature 
perturbations is regularized. It remains to study the coupling of $ to isocurva- 
ture perturbations case by case. We observe that this trick of a fictions total 
sound speed has been used in [13] in order to evolve $ and ignoring isocurvature 
perturbations. This of course is equivalent to study only half of the system and 
neglect the possibile interaction between isocurvature and adiabatic mode. 

The metric fluctuations induced by the fluid which drives the contraction $_|_ 
go through the bounce and seed the constant (growing) mode in the expansion. 
The metric fluctuations induced by the exotic fluid 4>_ are important only close 
to the bounce and decay with time in the expanding phase. However, we note 
that the latter have a spectrum which is steeper in the infrared than the constant 
mode. This means that the largest piece of $ in the contracting phase matches 
to a decaying mode in the expanding phase [21]. The picture which emerges 
seems the following: in presence of a second entity which is responsible for the 
bounce, both the growing and decaying mode of metric fluctuations seeded by 
the dominant fluid which drives the contraction match to the growing mode of 
$ in the expanding phase. 

In the expanding phase, after the bounce, gravitational waves have a scale 
invariant spectrum for a dust contraction as originally discovered by Starobinsky 
[17]. This is also what we find here (see also [16, 18]). In this two field model, 
after the bounce, scalar fluctuations with a (nearly, slightly red-tilted) scale 
invariant spectrum are generated by a very slow contraction, corresponding to 
a contracting fluid with an ultra-stiff equation of state {w+ >> 1). In this 
picture, the scale invariant spectrum for sub-Hubble metric fluctuations <I> is 
adiabatically transferred to super-Hubble wavelengths by the slow contraction, 
but remains imprinted in the subsequent expanding phase because of the bounce 
induced by a second entity. The amplitude of in agreement with observations 
could be obtained by tuning the parameters of the model. In the same way 
hopefully, should be kept in the linear regime during the bounce. This 
possibility could be of interest for the Ekpyrotic scenario [19], without the need 
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of a singularity. The model obtained would be similar to the idea of hybrid 
inflation in the context of inflationary models. As in hybrid inflation [20] one 

field drives inflation and the other terminates the accelerating stage, in this case 
one component would drive the contraction and the other would be responsible 
for the bounce. 

It is possible that the above conclusions arc due to the particular relation 
among the equations of state of the two fluids or to the time simmetry of the 
model, but this seems unlikely. The absence of intrinsic non adiabatic pressure 
perturbations (i. e. the assumption of two perfect fluids) may be a more im- 
portant issue. The danger could be a coupling (it may be just of gravitational 
origin) of to In such a case a scale invariance of metric fluctuations 
in the expanding phase may be very difficult to obtain. In any case this toy 
model shows how the physics at the bounce can alter the predictions about 
cosmological perturbations also in the PBB scenario [1]. 

It remains to connect this study with a prescription for matching perturba- 
tion through a non-singular bounce in a multifield case. 



Note Added: While this paper was almost completed, two papers which 
studied metric perturbations during a bounce appeared [22, 23]. Gasperini, 
Giovannini and Veneziano [22] studied a single field non singular bounce induced 
by a non local potential and found that the growing mode of $ in the contracting 
phase (not supported by a growing mode of C) does not match to the growing 
of $ in the expanding phase [24, 21]. Previous studies with a bounce induced 
by loop corrections (with or without a local potential) [25, 26] agree with [22]. 
This result could be related to the single field assumption. ToUey, Turok and 
Steinhardt [23] studied cosmological perturbations through a singular bounce. 
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